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1 Introduction 



One of the key issues of present high energy physics is to check whether 
the Higgs field(s) account for the mechanism of the mass generation. In the 
standard model (SM) of the high energy physics, only one physical Higgs 
boson exists after symmetry spontaneously broken, while five Higgs bosons 
in the minimal supersymmetric standard model (MSSM). 

Until now, the precise high energy experiments have tested the gauge 
bosons self interactions and the gauge bosons interactions with fermions (ex- 
cept top sector) of the SM with very good accuracy, while the Higgs sector 
is totally untouched [1]. In the SM and MSSM, the couplings among Higgs 
bosons and top quarks play a very unique role because which are not sup- 
pressed by ^ SL . Very likely, these couplings are firstly detected at experi- 
ments. Besides the aspect of the study on the properties of the Higgs bosons, 
there is another strong motivation to examine the H — t — t coupling: search- 
ing new physics beyond the SM. It is commonly thought that new physics 
might appear in the top sector for its huge mass. 

Supersymmetric (SUSY) QCD interactions by exchanging gluinos and 
squarks might impact greatly on the Higgs-top quarks vertexes. If the Higgs 
bosons are heavy enough, the SUSY QCD effects can be studied in the de- 
cay channel H — > tt (H is general notation for Higgs boson) [2]. Otherwise, 
production channels offer alternative methods and should be carefully exam- 
ined. At hadron colliders, Higgs boson can be associated production with top 
quark pair through gg and/or qq sub-processes [3]. Besides hadron colliders, 
as a supplementary machine, the linear colliders (LC) will provide ideal place 
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to study such effects. Even in the seventies, the process e + e~ — > Hit in the 
SM at LC has been proposed, which can be used as finding Higgs boson 
and measuring the couplings of gmt [4]. The calculations of the important 
QCD effects for the process are not accomplished until the end of the last 
century [5]. In the MSSM, the neutral Higgs boson associated production 
with top quark pair at LC including the QCD effect are available soon after 
[6]. However, the SUSY QCD effects on the process e + e~ — > tiH are not 
touched, which is the subject of this paper. In this paper, we focus on the 
lightest neutral Higgs boson h° associated production with top quark pair in 
the MSSM. 

2 Basics and technical set-ups 

The diagrams including the SUSY QCD corrections are shown in Fig. 1-3. 
These diagrams and Feynman amplitudes are generated by FeynArts [7]. 
Fig. 1 and 2 are the Feynman diagrams at tree and one-loop level, respec- 
tively. In Fig. 3, we draw the diagrams containing counter-terms. In this 
paper, we renormalize the quark wave function and mass in the on-shell (OS) 
renormaliztion scheme. Therefore in Fig. 3, we do not show the diagrams 
with counter-terms at external legs because in the OS renormaliztion scheme, 
such kind of diagrams are completely eliminated by the corresponding virtual 
diagrams. 

The analytical results for Feynman amplitudes and phase space integra- 
tions are tedious but easily obtained through standard procedure, which are 
useless to almost all readers. Therefore, we decide to omit the explicit analyt- 
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ical results here 2 . For the calculations of the amplitudes and the Monte Carlo 
phase space integrations, we rely heavily on the popular packages FormCalc 
and LoopTools [8] 3 . Before we shift to numerical discussions, it should be 
emphasized that the entire calculations and the cancellation of ultraviolet 
divergences have been checked carefully. 

3 Numerical discussions and summary 

In the MSSM, there are dozens of parameters and experiments constrains, 
for examples Ap and direct mass limits for gluino, squarks and Higgs bosons. 
We rely on the built-in programs of FormCalc to calculate the SUSY par- 
ticles mass spectrum (we use two-loop results for Higgs mass calculations) 
and SUSY parameters as well as set constraints on known limits from exper- 
iments. Specifically, we set 

Ap < 3 x KT 3 (1) 

m-g > 175 GeV, m~ t > 80 GeV, m~ b > 70 GeV 
m~ q > 150C7eU, m h o > 91GeV. (2) 

For simplicity we have assumed the SUSY parameters Mq = M v = M D = 
M s and A t = A b . 

As for the first attempt, we would like to see the SUSY QCD effects for 

the non-mixing case, i.e. the non-diagonal matrix elements Mlr = A t — 
2 Programs are available on request. 

3 A technical detail, LoopTools is not applicable for some Do with certain parameters 
that all (mf + m j ~Pij) 2 < ^mfm?, details of the definition and convention for D can be 
found in manual of LoopTools or in Ref [9]. For this exceptional case, the new program 
has been added according to Ref [9]. 
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///tan/3 for stop mass matrix are zero. Our numerical results show that 
the SUSY QCD effects are small for such case, and the relative corrections 
5 = ((XjVLO - vlo)/°'lo are at most about two percents, which are hard to 
detect at experiments. In the following we will focus on the mixing case. 

After scanning the whole parameters space, in order to demonstrate the 
numerical results, in Fig. 4-6 we show the cross sections and relative cor- 
rections 5 for a specific parameter set as a function of tan f3, center-of-mass 
energy y/s of e + e~, as well as the mass of gluino. From the figures, we can 
see that the next-to-leading order (NLO) SUSY QCD effects decrease the 
LO cross sections and the magnitude can reach 50%. For such kind of large 
corrections, it should be noted that main contributions come from diagrams 
(5)- (8) of Fig. 2, and the lighter stop mass is about 100 GeV in this case. 
Similar to the SM QCD effects, the large corrections occur in the vicinity of 
the production threshold (Fig. 5), where the cross sections become small. 
Furthermore, from Fig. 6, it is obviously that the corrections are at the 
peak around trig = 190 ~ 200 GeV, where the re-scattering enhancement 
t* — > g + ti — > t + h° appears at rrig + ~ m t + mho. 

In table 1, we show that cross sections at LO and NLO as well as the rela- 
tive corrections as a function of A t . The large SUSY QCD effects emerge only 
with positive large A t [or Mlr], which ensures that the large contributions 
from lighter stops don't vanish. 

From our parameters space scanning, in order to have the large SUSY 
QCD effects, the following requirements should be satisfied 

1. there must be large mass splitting between stops, moreover Mlr should 
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be large and positive, 

2. M s should not be large, 

3. gluino mass should not be heavy. 

As mentioned above, condition (1) prevents the cancellation between stops, 
and makes one stop mass light which is essential for large SUSY QCD effects, 
especially in the threshold regions. Conditions (2) and (3) ensure that SUSY 
effects do not decouple from the process. 

To summarize, the SUSY QCD corrections to the lightest neutral Higgs 
boson associated production with top quark pair are studied in the MSSM 
at LC. Our studies show that generally the corrections are very moderate 
(say 10%) and under control, except for some rescattering effects which lead 
to a breakdown of perturbation theory and require a more detailed study, 
resummation or something of this kind. Contrary to the QCD corrections 
arising from gluons which are about 50% for m h o = 120 GeV and y/s = 
500 GeV [near the production threshold], the SUSY QCD effects can be 
also as large as about -50% in this kinemical region for the favorable model 
parameter space under OS renormalization scheme. Such kind of large SUSY 
QCD effects might be the probe to determine the sign of M LR . We note that if 
nature does choose lighter stops, which can be copiously produced at hadron 
colliders and linear colliders. Comparison between direct productions and 
indirect virtual effects might give more information on the MSSM. 
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A t (GeV) 


1000 


950 


900 


800 


700 


500 


200 


myfi (GeV) 


113.4 


115.2 


116.1 


115.7 


113.5 


106.8 


97.6 


m h (GeV) 


99.4 


136.4 


165.3 


211.5 


249.4 


311.5 


386.4 


LO 


0.463 


0.411 


0.389 


0.400 


0.459 


0.674 


1.075 


NLO 


0.198 


0.327 


0.348 


0.387 


0.452 


0.670 


1.077 


S(%) 


-57.2 


-20.4 


-10.4 


-3.2 


-1.6 


-0.7 


0.2 


A t (GeV) 





-200 


-500 


-700 


-800 


-900 


-1000 


m h o (GeV) 


94.7 


96.6 


105.1 


112.0 


114.7 


116.2 


115.2 


m h (GeV) 


429.1 


397.5 


325.2 


266.3 


231.2 


189.8 


136.4 


LO 


1.226 


1.128 


0.741 


0.503 


0.425 


0.387 


0.411 


NLO 


1.229 


1.131 


0.741 


0.506 


0.432 


0.393 


0.410 


n%) 


0.3 


0.3 





0.5 


1.7 


1.6 


-0.2 



Table 1: Total cross sections [fb] at LO and NLO as well as the relative 
corrections as a function of A t with tan/3 = 4. Other parameters are the 
same with Fig. 4. 
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Figure 1: Born diagrams of e + e — > 
9 
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Figure 4: Cross sections (a), relative correction (b), Higgs mass (c) and 
lighter stop mass (d) as a function of tan/5 for ^/s = 0.5 TeV with A t — 1 
TeV, M s = 400 GeV, /i = -100 GeV, m A o = 300 GeV, m § = 200 GeV. The 
solid and dashed lines in (a) reDresent the cross sections at NLO and LO. 
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Figure 5: Cross sections (a) and relative correction (b) as a function of ^/s for 
tan/5 = 4, where m h o = 113.3 GeV and = 99.4 GeV. Other parameters 
and conventions are the same with Fig. 4. 




14 



Figure 6: Cross sections (a) and relative correction (b) as a function of 
rrig with tan/5 = 4, where m^o = 113.3 GeV and = 99.4 GeV. Other 
parameters and conventions are the same with Fig. 4. 



